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Terahertz  Radiation:  A  Non-contact  Tool  for  the 
Selective  Stimulation  of  Biological  Responses 

in  Human  Cells 

Ibtissam  Echchgadda,  Jessica  E.  Grundt,  Cesario  Z.  Cema,  Caleb  C.  Roth,  Jason  A.  Payne, 

Bennett  L.  Ibey,  and  Gerald  J.  Wilmink 


Abstract — Collective  motions  of  water  and  of  many  biological 
macromolecules  have  characteristic  time  scales  on  the  order  of  a 
picosecond.  As  a  result,  these  biomolecules  can  strongly  absorb 
terahertz  (THz)  radiation.  Due  to  this  absorption,  THz  radiation 
can  exert  a  diverse  range  of  effects  on  biological  structures.  For 
example,  THz  radiation  has  been  shown  to  impact  the  structure, 
functional  activity,  and  dynamics  of  macromolecules  such  as  DNA 
and  proteins.  THz  radiation  can  affect  several  gene  expression 
pathways  and,  consequently,  can  alter  various  biochemical  and 
physiological  processes  in  cells.  Indeed,  THz  radiation  has  been 
shown  to  influence  the  expression  of  several  genes  within  different 
cell  types.  However,  a  complete  view  of  the  global  transcriptional 
responses  and  the  intracellular  canonical  pathways  specifically 
triggered  by  THz  radiation  has  not  been  elucidated.  In  this  study, 
we  performed  a  global  profiling  of  transcripts  in  human  cells  ex¬ 
posed  to  2.52  THz  radiation  and  compared  the  exposure  responses 
to  a  thermally-matched  bulk-heating  (BH)  protocol.  Our  results 
show  that  both  THz  radiation  and  BH  induce  a  significant  change 
in  the  expression  of  numerous  mRNAs  and  microRNAs.  The  data 
also  show  that  THz  radiation  triggers  specific  intracellular  canon¬ 
ical  pathways  that  are  not  affected  in  the  BH-exposed  cells.  This 
study  implies  that  THz  radiation  may  be  a  useful,  non-contact  tool 
for  the  selective  control  of  specific  genes  and  cellular  processes. 

Index  Terms — Cell  exposure,  cellular  response,  global  gene 
expression,  miRNAs,  mRNAs,  signaling  pathways,  thermal  effect, 
THz  bioeffect. 


I.  INTRODUCTION 

TERAHERTz  (THz)  radiation  strongly  interacts  with 
biomolecular  systems  whose  time  scales  are  on  the  order 
of  a  picosecond  (ps).  Recent  experimental  and  molecular  dy¬ 
namics  simulations  indicate  that  bulk  water,  interfacial 
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water  (water  molecules  positioned  at  the  interface  of  macro¬ 
molecules),  and  numerous  macromolecules  (e.g.,  carbohy¬ 
drates,  proteins,  and  nucleic  acids)  strongly  absorb  electromag¬ 
netic  radiation  with  frequencies  ranging  from  0.1  to  10  THz, 
where  1  THz  is  equivalent  to  1  ps-1. 

Water  molecules  possess  potent  absorption  of  THz  radiation 
due  to  their  strong  dipole  moment  and  their  high  intermolec- 
ular  hydrogen-bond  (HB)  rearrangements  (i.e.,  intermolecular 
stretching  and  bending  vibrational  modes)  [1]— [5].  Bulk  water 
collective  HB  networks  have  been  shown  to  exhibit  intermolec¬ 
ular  stretching  and  bending  vibrational  modes  at  5  and  1 .5  THz, 
respectively  [1],  [2],  [5].  Interfacial,  also  termed  solvation, 
water,  which  plays  a  central  role  not  only  in  stabilizing  the 
structure  of  macromolecules  but  also  in  mediating  their  specific 
function  [6]— [1 1],  exhibits  dynamics  and  vibrational  motion 
that  differ  considerably  from  that  of  bulk  water.  Specifically, 
bulk  water  consists  of  four  HBs  that  rapidly  rearrange  on  the 
ps  timescale,  whereas  interfacial  water  predominantly  forms 
two  HBs  along  the  axis  of  a  biomacromolecular  solute,  which 
makes  HBs  in  interfacial  water  move  slower  and/or  survive 
longer  than  those  in  bulk  water  [12].  Recent  THz  spectroscopy 
and  molecular  dynamics  simulations  indicate  that  interfacial 
water  creates  two  distinct  hydration  shells,  which  extend  up 
to  20  A  from  the  surface  of  the  biomolecule  [12]— [14].  These 
studies  also  indicate  that  the  second  solvation  shell  exhibits  a 
prominent  absorption  peak  between  2.4  and  2.5  THz  [12]  and 
that,  over  the  2.2-2.9-THz  frequency  range,  interfacial  water 
exhibits  a  higher  absorption  coefficient  (//a  «  400  cm-1) 

than  bulk  water  [12].  Molecular  dynamics  simulations  confirm 
that  these  differences  are  a  result  of  the  difference  in  structure 
between  interfacial  and  bulk  water.  In  addition  to  the  interfacial 
water  molecules  that  solvate  them,  macromolecules  also  exhibit 
high  absorption  coefficients  at  THz  frequencies  largely  because 
of  their  strong  low  frequency  collective  modes  (i.e.,  skeletal 
and  breathing  modes),  which  occur  on  the  ps  time  scale  [12], 
[13],  [15]-[22]. 

Due  to  the  strong  absorption  of  THz  radiation  by  water  and 
macromolecules,  computational  models  predict  that  exposure 
to  high-power  THz  radiation  will  cause  an  appreciable  amount 
of  energy  buildup  within  the  exposed  cells.  This  localized 
energy,  while  too  low  (in  the  meV  range)  to  cause  ionizing 
damage,  is  predicted  to  result  in  significant  local  temperature 
rise  leading  to  significant  effects  at  the  cellular,  organelle,  and 
macromolecular  levels  [23]-[28].  Recent  studies  investigating 
the  biological  effects  of  THz  radiation  have  shown  that  this  radi- 
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ation  interacts  with  the  plasma  membrane  [24],  [25],  impacting 
its  organization,  function,  and  permeability  [24]-[28].  Several 
in  vitro  studies  performed  on  biological  macromolecules  have 
shown  that  THz  radiation  can  also  directly  affect  the  structure, 
functional  activity  and  dynamics  of  DNA  and  proteins  [26], 
[29]— [3 1  ]. 

THz  radiation  can  also  significantly  change  the  normal 
process  of  gene  expression  and,  consequently,  can  influence  the 
cellular  phenotypic  properties  (e.g.,  morphology,  biochemical 
or  physiological  characteristics).  In  fact,  a  number  of  reports 
have  showed  that  both  exposure  to  a  continuous-wave  (CW) 
THz  radiation  and  THz  radiation  from  a  pulsed  broadband 
source  can  cause  gene  transcriptional  alterations  in  skin  tissues 
[32],  [33]  and  various  cell  types  [25],  [34]— [39].  The  recent 
experiments  with  broadband  THz  pulses  were  conducted  under 
controlled  thermal  conditions  and  the  temperature-increases 
were  measured  to  be  not  significant.  Therefore,  these  studies 
suggested  that  THz  radiation  can  interact  with  biomolecules 
and  cells  through  non-thermal  mechanisms  (resonance-type 
linear  and  nonlinear  interactions).  They  also  argued  that  the 
shift  in  gene  transcriptional  activity  by  THz  radiation  is  in  part 
due  to  a  change  in  the  molecular  function  of  macromolecules. 
For  example,  THz  irradiation  may  have  directly  caused  an 
alteration  in  the  breathing  vibrational  dynamic  of  DNA  and 
proteins  leading  to  a  change  in  their  shapes  or  conformational 
states  [40]— [45],  and  thereby  affecting  chromatin  organi¬ 
zation,  hindering  or  inducing  protein-DNA  binding,  and/or 
protein-protein  interactions  between  transcription  factors  and 
coregulators. 

Several  of  the  groups  employing  high-power  CW  THz  radia¬ 
tion,  on  the  other  hand,  believe  that  given  the  high  absorption  of 
this  radiation  by  water  in  biological  tissue,  the  THz-induced  ef¬ 
fects  on  gene  expression  can  be  mediated  by  heat  (hyperthermic 
effects),  which  results  from  the  temperature  buildup  during  ex¬ 
posure  [26],  [28],  [3 1].  Therefore,  it  is  proposed  that  the  temper¬ 
ature  rises  can  cause  deregulation  due  to  unfolding,  entangle¬ 
ment,  and  aggregation  of  intracellular  and  membrane  proteins 
[46]— [53].  In  this  light,  the  question  of  whether  cells  exposed 
to  intense  CW  THz  radiation  exhibit  similar  cellular  effects  as 
those  exposed  to  a  matched  conventional  thermal  bulk-heating 
(BH)  is  important. 

In  this  study,  we  aimed  to  identify  the  mRNAs,  microRNAs 
(miRNAs),  and  the  intracellular  cascades  that  are  specifically 
triggered  by  high-power  CW  THz  radiation.  Previously  con¬ 
ducted  experiments  have  demonstrated  that  a  CW  single-fre¬ 
quency,  2.52  THz  radiation  can  induce  a  change  in  the  expres¬ 
sion  of  several  genes  [31],  [34]— [36].  However,  these  studies 
fail  to  provide  a  comprehensive  view  of  the  global  genomic 
and  cellular  responses  triggered  by  THz  irradiation.  Addition¬ 
ally,  they  fail  to  characterize  THz  radiation  transcriptional  effect 
on  the  non-coding  miRNAs.  This  characterization  is  significant 
given  that  miRNAs  are  considered  master  regulators  of  gene  ex¬ 
pression  and  are  critical  players  in  many  fundamental  cellular 
processes  including  development,  proliferation,  differentiation, 
death,  and  metabolism  [54]— [58]. 

We  conducted  experiments  where  human  Jurkat  T  lympho¬ 
cytes  were  exposed  to  a  CW,  single-frequency,  2.52  THz  or 
a  thermally  matched  BH  protocol.  Global  expression  profiling 


Fig.  1.  Experimental  setup  for  exposure  of  cells  to  2.52  THz  radiation, 
(a)  Schematic  representation  of  our  custom  THz  exposure  system,  (b)  Magni¬ 
fication  of  THz  transmission  and  delivery  optics,  (c)  Magnification  of  a  single 
well  with  suspension  cells  pooled  at  bottom. 


of  transcripts  and  cellular  responses  were  evaluated  using  mi¬ 
croarray  and  bioinformatics  analysis  tools.  Overall,  our  results 
demonstrated  that  THz  radiation  induced  a  significant  change  in 
the  expression  of  specific  mRNAs  and  miRNAs.  Furthermore, 
the  data  also  show  that  THz  radiation  triggered  specific  intracel¬ 
lular  metabolic  and  signaling  pathways  that  were  not  identified 
in  the  BH-exposed  cells.  Overall,  these  data  provide  valuable 
new  insights  that  give  a  much  clearer  picture  of  the  genes  and 
intracellular  pathways  that  are  specifically  triggered  in  human 
cells  exposed  to  THz  radiation.  Additionally,  the  data  suggest 
that  THz  radiation  may  be  a  useful,  non-contact  tool  for  the  se¬ 
lective  control  of  specific  genes  and  cellular  processes. 

II.  Methodology 

A.  Custom-Designed  THz  Exposure  Enclosure 

To  conduct  temperature-controlled  THz  radiation  exposures, 
we  used  our  custom-designed  exposure  system  [31],  [34], 
which  consists  of  a  temperature  controller  (set  to  37  °C), 
nitrogen  gas  purging  system,  operator  control  area  equipped 
with  polyvinyl  chloride  (PVC)  gloves,  input  ports  for  both 
short  and  long  THz  radiation  cavities,  optical  platform,  and 
a  sapphire  window  for  infrared  (IR)  thermographic  measure¬ 
ments.  Fig.  l(a)-(c)  represents  an  illustration  of  our  custom 
THz  exposure  setup.  To  ensure  the  chamber  was  exhibiting 
consistent  heat  distribution,  Omega  thermocouples  were  used 
to  monitor  the  chamber’s  thermal  stability  at  various  spatial 
positions.  Temperature  readings  were  collected  from  the  bottom 
to  the  top  of  the  chamber  using  2  cm  increments  (z-axis).  The 
data  showed  that  the  exposure  chamber  maintained  consistent 
optimal  temperatures  (<  0.5  °C)  at  all  positions. 

B.  THz  Beam  Manipulation  and  Delivery  Optics 

For  all  exposures,  we  used  an  optically  pumped  molecular  gas 
THz  radiation  laser  source  (SIFIR-50  OPTL,  Coherent-DEOS, 
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Fig.  2.  (a)  Plot  of  power  intensity  as  function  of  frequency  at  THz.  (b)  Sample 
representative  image  of  THz  intensity  beam  profile,  (c)  FDTD  model  results  of 
the  THz  exposure  of  cellular  medium  within  the  384-well  Cyclo-Olefin  polymer 
plate  geometry. 


Santa  Clara,  CA,  USA)  [31],  [34].  Three  optical  elements  were 
used  to  control  and  deliver  the  THz  beam  from  the  source  to 
the  exposure  plate.  The  first  element  is  a  microprocessor-con¬ 
trolled,  motor-driven  shutter  (Sutter  Instr.  Co,  Novato,  CA). 
This  shutter  has  a  25-mm  diameter  aperture  and  a  38-ms  open/ 
close  time.  To  eliminate  appreciable  cutoff  due  to  THz  beam  di¬ 
vergence,  the  shutter  was  positioned  7  cm  from  the  opening  of 
the  short-wave  cavity.  Once  delivered  through  the  shutter  aper¬ 
ture,  a  flat  gold  mirror  (r  =  25  mm)  and  a  parabolic  silver 
mirror  were  used  to  collect  and  focus  the  incident  THz  beam 
on  the  well  plate. 


C.  THz  Detectors  and  Beam  Diagnostics 

Prior  to  conducting  cellular  exposures,  comprehensive  laser 
beam  diagnostics  were  performed  to  ensure  the  THz  source  was 
delivering  desirable  THz  radiation.  THz  output  power  was  mea¬ 
sured  with  an  Astral  Vector  Series  H410  calorimeter  (Scien- 
tech,  Boulder,  CO).  Spatial  beam  profiles  were  evaluated  using 
a  Spiricon  Pyrocam  III  detector  array  (Spiricon  Inc,  Logan,  UT, 
USA).  As  seen  in  a  sample  image  in  Fig.  2(a),  we  found  that  the 
THz  radiation  source  provided  an  irradiance  of  636  mWcm2  at  a 
center  frequency  of  2.52  THz  and  a  full-width  at  half-maximum 
(FWHM)  of  0.04  THz.  Using  a  pyroelectric  array  camera,  we 
found  that  the  intensity  beam  profile  exhibits  Gaussian  distribu¬ 
tion  with  a  diameter  of  3.4  mm,  as  shown  in  Fig.  2(b). 


D.  Dosimetry:  Temperature  History  Before,  During,  and  After 
THz  Exposure 

Jurkat  cells  (TIB- 152,  ATCC,  Manassas,  VA),  a  human  T 
lymphocyte  cell  line,  were  maintained  in  Advanced  RPMI 1 640 
media  (Life  Technologies,  Grand  Island,  NY)  supplemented 
with  10%  fetal  bovine  serum,  10  mM  Hepes,  100  units/ml 
penicillin,  and  100  /ig/ml  streptomycin  in  a  humidified,  5% 
CO2  incubator  at  37  °C.  For  all  THz  radiation  exposure  ex¬ 
periments,  cells  were  plated  in  384-well  Cyclo-Olefin  tissue 
culture  plates  (Aurora  Biotechnologies,  Carlsbad,  CA)  at  a 
seeding  density  of  ~  1.6  x  105  cells  per  well  (140  /il  total 
volume  of  cells  and  medium).  Cells  were  incubated  overnight 
and  exposed  the  next  day,  as  previously  described  [31].  The 
Aurora  384-well  Cyclo-Olefin  polymer  tissue  culture  plates  are 
equipped  with  an  evaporation  control  system  that  minimizes 
gas  exchange.  In  addition,  for  added  precaution,  prior  to  THz 
exposures  we  further  sealed  the  top  of  each  well  plate  with  an 
optical  adhesive  film  to  prevent  any  suspected  evaporation  of 
media.  Cells  were  exposed  to  THz  radiation  (y  =  2.52  THz, 

=  636  mWcm— 2)  from  below  for  40  min.  The  main  goal 
for  our  initial  set  of  experiments  was  to  determine  the  thermal 
history  for  Jurkat  cells  exposed  to  THz  radiation.  To  this  end, 
we  used  computational  and  empirical  dosimetric  tools  [31]. 

For  the  computational  dosimetry,  we  employed  finite-dif¬ 
ference  time-domain  (FDTD)  modeling  techniques  to  predict 
the  incident  fields  at  the  cells  while  accounting  for  the  well 
plate  geometry  that  replicate  the  empirical  exposure  conditions. 
Fig.  2(c)  represents  our  modeling  of  the  THz  exposure  of  cel¬ 
lular  medium  within  the  384-well  Cyclo-Olefin  polymer  plate 
geometry.  Within  this  simulation,  the  well  plate  geometry  was 
modeled  by  cellular  medium  (relative  permittivity  =  4.4  and 
conductivity  =  169.8  Siemens  m-1)  surrounded  by  lossless 

polystyrene  walls  (relative  permittivity  =  2.6)  mounted  on  top 
of  a  lossless  Cyclo-Olefin  polymer  film  (relative  permittivity 
=  2.3).  The  results  of  this  simulation  indicated  an  incident 
field  intensity  of  532  mW/cm2  at  the  location  of  the  cell 
layer,  corresponding  to  approximately  1 6%  reduction  from  the 
636  mW/cm— 2  measured  with  the  calorimeter.  These  simulated 
field  values  were  then  inputted  into  a  finite-difference  thermal 
solver  to  predict  the  change  in  temperature  over  time  due  to  the 
THz  exposure,  as  described  before  [31],  [34]. 

For  empirical  temperature  measurements,  we  used  IR  ther¬ 
mography  and  thermocouples.  The  IR  camera  was  placed 
1 5  cm  from  the  media  surface,  and  temperature  measurements 
were  acquired  at  800  frames/s  (SC6000,  FLIR  Systems  Inc, 
Wilsonville,  OR,  USA).  The  thermocouples  had  the  following 
specifications:  T-type,  Copper  material,  0.005”  diameter,  Teflon 
insulation,  and  a  response  time  of  40  ms  (Omega  Engineering, 
Inc.  Stamford,  CT,  USA).  Temperature  data  were  collected 
with  thermocouples  before,  during,  and  after  THz  exposure. 
For  comparison,  the  temperature  changes  in  sham  (unexposed) 
control  were  also  recorded.  The  sham  control  was  placed 
near  the  THz-exposed  cells  in  the  custom-designed  37  °C 
preconditioned  THz  exposure  chamber.  Sample  representative 
time-lapse  IR  images  are  provided  in  Fig.  3(a)  and  a  compara¬ 
tive  graph  that  depicts  the  thermal  data  for  both  empirical  and 
computational  methods  is  shown  in  Fig.  3(b). 


ECHCHGADDA  et  al.\  THZ  RADIATION:  NONCONTACT  TOOL  FOR  THE  SELECTIVE  STIMULATION  OF  BIOLOGICAL  RESPONSES  IN  HUMAN  CELLS 


57 


0  10  20  30  40  50 

THz  exposure  duration  (minutes) 


Fig.  3.  Temperature  history  of  cells  exposed  to  2.52-THz  radiation,  (a)  Time- 
lapse  IR  images  of  temperature  profiles  for  the  surface  of  media  in  a  well  before, 
during,  and  after  exposure,  (b)  Computational  and  quantitative  temperature-time 
data  collected  for  a  well  exposed  to  THz  radiation  using  finite-difference  time- 
domain  (FDTD)  modeling  techniques  and  thermocouples,  respectively. 


Overall,  although  our  computational  FDTD  model  shows  that 
the  temperature  of  the  well  rises  sharply  and  more  quickly  com¬ 
pared  to  the  empirical  data  (thermocouples),  they  both  flatten 
out  at  about  similar  level  for  the  remainder  of  the  exposure.  Both 
techniques  were  in  agreement,  showing  an  average  temperature 
rise  of  ~  6  °C  during  THz-cell  exposure.  This  change  in  tem¬ 
perature  (AT  °C)  value  was  used  to  set  the  thermally  matched 
conditions  of  the  BH.  Thus,  the  BH-exposed  cells,  seeded  as 
the  THz  exposed  cells,  in  384-well  Cyclo-Olefin  polymer  tissue 
culture  plates  with  sealed  tops,  were  placed  in  a  controlled  water 
bath  fixed  at  44  °C  for  40  min.  This  temperature  was  selected  to 
match  the  final  temperature  observed  in  the  THz  radiation  ex¬ 
posed  cells. 

Additionally,  we  conducted  MTT  (3-(4,5-Dimethylthiazol-2- 
yl)-2,5-diphenyltetrazolium  bromide)  assays  to  determine  the 
viability  of  cells  exposed  to  THz  radiation  or  BH.  Cell  viability 
was  evaluated  using  MTT  assays,  as  per  manufacturer's  instruc¬ 
tions  (ATCC,  Manassas,  VA,  USA).  Our  experiment  settings  for 
cell  exposure  to  THz  radiation  or  BH  resulted  in  a  similar  cel¬ 
lular  viability  measure,  which  was  about  70%  compared  to  sham 
cells. 


E.  RNA  Isolation  and  Normalization 

mRNAs  and  miRNAs  were  harvested  from  the  sham,  THz 
radiation-  and  BH-exposed  cells  using  the  RNeasy  Mini-Kit 
and  the  miRNeasy  mini  kit  (Qiagen,  Valencia,  CA,  USA), 
respectively,  as  per  the  manufacturer's  instructions.  RNA 
extraction  was  performed  at  4  h  post- expo  sure.  The  RNAs' 
concentrations  were  assessed  on  a  NanoDrop  Spectropho¬ 
tometer  (NanoDrop  Technologies,  Wilmington,  DE,  USA), 
and  the  quality  was  measured  on  a  2100  Bioanalyzer  (Agilent 
Technologies,  Santa  Clara,  CA).  RNA  samples  with  a  RNA 
Integrity  Number  greater  than  9.5  were  subjected  to  mRNA 
and  miRNA  microarray  analysis. 


F  mRNA  Microarrays  and  PCR 

To  determine  the  effect  of  THz  radiation  and  BH  exposures 
on  mRNA  expression,  we  conducted  microarray  analysis  for 
each  exposure  condition  (sham,  THz  radiation,  and  BH)  in 
triplicate  using  the  Affymetrix  GeneChip  Human  Genome 
U133  (HG-U133  plus  2.0  Array  that  contains  54  675  probe  sets, 
Asuragen  Services,  Austin,  TX,  USA).  Briefly,  two  micrograms 
of  RNA  were  used  for  preparation  of  biotin-labeled  targets 
(cRNA)  using  MessageAmp-based  protocols  (Ambion,  Inc.). 
Labeled  cRNA  was  fragmented  and  used  for  array  hybridiza¬ 
tion  and  washing.  The  cRNA  was  mixed  with  a  hybridization 
cocktail,  heated  to  99  °C  for  5  min,  and  then  incubated  at 
45  °C  for  5  min.  Hybridization  arrays  were  conducted  for 
16  h  in  an  Affymetrix  Model  640  hybridization  oven  (45  °C, 
60  rpm).  Arrays  were  washed  and  stained  on  a  FS450  Fluidics 
station  and  scanned  on  a  GeneChip  Scanner  3000  7G.  Image 
signal  data,  detection  calls,  and  annotations  were  generated  for 
every  gene  using  the  Affymetrix  Statistical  Algorithm  MAS 
5.0  (GeneChip  Operating  Software  vl.3).  A  log2  transforma¬ 
tion  was  conducted  and  a  Student's  t-test  was  performed  for 
comparison  of  the  two  groups  (sham  and  exposed).  Multiple 
testing  corrections  using  Benjamini  and  Hochberg  methods 
were  used  to  determine  the  false  discovery  rate,  and  genes  with 
statistically  significant  results  were  identified  using  Bonferroni 
correction  procedures  (—  log  10  pcutoff  >  6.04)  [59].  For 
interpretation  of  the  results,  the  Ingenuity  Pathways  Analysis 
tool  (IPA,  version  8.7,  Ingenuity  Systems  Inc.,  Redwood  City, 
CA,  USA1)  was  used.  IPA  is  a  web-based  software  applica¬ 
tion  that  enables  filtering  and  dataset  comparisons  to  identify 
biological  mechanisms  pathways  and  functions  most  relevant 
to  experimental  datasets  or  differentially  expressed  genes. 
The  cutoff  criteria  for  IPA  analysis  were  an  absolute  value  of 
log2  ratio  >  2  and  a  p  —  value  <  0.05.  Other  web-based 
resources  were  used  to  complement  the  interpretation  of  the 
data.  These  included  the  GeneCards  Human  Gene  Database 
(http://www.genecards.org/)  and  The  HUGO  Gene  Nomencla¬ 
ture  Committee  (HGNC,  http://www.genenames.org/).  A  few 
selected  target  mRNAs  identified  in  the  microarray  analysis 
were  validated  using  quantitative  real-time  polymerase  chain 
reactions  (qRT-PCR).  Runs  were  performed  on  a  StepOnePlus 
system  using  TaqMan  RNA-to-CT  1-Step  Kits  and  TaqMan 
Assays  (Applied  Biosystems  by  Life  Technologies).  Calibrator 
RNA  was  used  as  a  control.  PCR  was  conducted  using  a 
three-program  LightCycler  protocol. 

G.  miRNAs  Microarray  Analysis  and  PCR 

For  the  miRNA  microarrays,  an  Affymetrix  GeneChip  (Dis- 
covArrayTM,  Asuragen  Services)  was  used,  which  included 
probes  of  all  miRNAs  from  the  Sanger  miRBase  and  more 
than  12,000  predicted  miRNAs  [60]-[66].  For  each  probe,  a 
two-sample  t-test  for  empirical  Bayes  variance  estimates  was 
used.  Probes  with  p  —  values  <  0.05  were  determined  to  be 
statistically  significant.  Normalized  log2  transformed  intensity 
values  were  analyzed  using  JMP  Genomics  3 .2  (S  AS  Institute, 
Cary,  NC).  IPA  was  used  for  further  refinement  and  interpreta¬ 
tion  of  miRNA  microarray  data.  The  analysis  cutoffs  were  set 

1  [Online].  Available:  http://www.ingenuity.com 


58 


IEEE  TRANSACTIONS  ON  TERAHERTZ  SCIENCE  AND  TECHNOLOGY,  VOL.  6,  NO.  1,  JANUARY  2016 


(a) 


£ 

jy- 

fi  , 

E'l  ’’ 

*  £  U  - 
S  —  M 

I  J 


III/ 


Bit 


£  ’ 


-toSw 

(p-vutur) 


ft  |  *  i*  ,IHl|  * 

*  A  *  >  .  .si*  i ■  -■■  r i , 


-  l«E,u 
(P-miIlM 


(C) 


□  up-regulated 


Fig.  4.  Differential  expression  profiles  of  several  mRNAs  in  THz  and  BH  ex¬ 
posed  cells,  (a)  Volcano  plots  of  the  gene  expression  profiles,  where  the  mag¬ 
nitude  of  differential  expression  (log2  fold-change)  is  plotted  versus  the  level 
of  statistical  significance  (—  log10  p-value).  Up-regulated  and  down-regulated 
genes  with  absolute  log2  ratios  >  2  and  p  —  values  <  0.05  are  depicted 
with  green  and  red  triangle  symbols,  respectively.  Genes  denoted  with  light  gray 
symbols  were  not  significantly  different,  (b)  qRT-PCR  validation  of  the  change 
in  expression  of  selected  mRNAs.  (c)  Total  number  of  up-  and  down-regulated 
genes  for  THz  and  BH  exposures.  Using  IPA,  the  data  show  that  531  and  672 
genes  were  differentially  expressed  in  THz  and  BH,  respectively. 


atp  —  values  <  0.05  and  an  absolute  value  of  log2  ratio  >  0.59 
(equivalent  to  1.5  fold  change  in  expression). 

III.  Results 

A.  Global  mRNAs  Expression  Profiles  for  Cells  Exposed  to 
THz  Radiation  or  BH 

To  get  a  comprehensive  view  of  the  transcriptional  response 
to  THz  radiation,  Jurkat  cells  were  exposed  to  2.25  THz  radi¬ 
ation  or  BH  and  microarray  gene  chips  were  used  to  quantify 
the  expression  levels  of  mRNAs  in  each  exposure  condition 
versus  control  sham  (unexposed  cells).  Since  the  microarrays 
provide  data  for  nearly  60,000  mRNA  targets,  volcano  plots 
were  created  to  best  visualize  the  global  distribution  of  the 
mRNAs  expression  profiles  for  THz  radiation  and  BH,  as 
presented  in  Fig.  4(a).  In  these  plots,  the  level  of  statistical 
significance  (p-value)  for  each  mRNA  is  plotted  against  the 
fold  change  of  expression  magnitude  relative  to  sham.  mRNAs 
with  the  highest  level  of  statistical  significance  appear  on  the 
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Fig.  5.  Identification  of  unique  and  shared  mRNAs.  (a)  Venn  diagram  indi¬ 
cating  the  number  of  unique  and  common  differentially  expressed  mRNAs  in 
each  exposed  group,  (b).  Gene  family  of  the  unique  and  shared  genes.  The  pie 
chart  sectors  represent  the  number  of  genes  associated  with  each  family  as  in¬ 
dicated  in  the  legend.  The  protruded  sectors  represent  the  four  gene  families 
(cytokine,  growth  factor,  peptidase,  and  translation  regulator)  that  did  not  dis¬ 
play  shared  genes  between  THz  and  BH  exposures. 


right  side  of  the  plot.  Significantly  up-regulated  mRNAs  appear 
towards  the  top  of  the  plot  (denoted  with  green  triangles), 
whereas  significantly  down-regulated  mRNAs  appear  towards 
the  bottom  (denoted  with  red  triangles).  The  gray  color  repre¬ 
sents  the  mRNAs  that  were  considered  unchanged.  Overall,  the 
global  changes  in  mRNA  expression  profiles  were  similar  for 
each  exposure  condition  since  1162  and  1200  mRNA  probes 
were  differentially  expressed  in  the  THz  radiation  and  BH 
exposure  conditions,  respectively.  Fig.  4(b)  shows  a  validation 
of  the  changes  in  mRNAs  expression  of  few  mRNAs  probes 
in  the  array  using  qRT-PCR.  These  data  show  that  the  mRNAs 
expression  levels  were  in  agreement  with  those  obtained  in  the 
microarray  study. 

IPA  software  was  used  to  refine  the  microarray  data.  Using  a 
p  —  value  <0.05  and  an  absolute  log2  ratio  >  2  (equivalent 
to  four- fold  change)  cutoffs,  531  and  672  distinct  mRNAs  were 
identified  as  differentially  expressed  in  THz  radiation  and  BH 
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Fig.  6.  Cellular  and  molecular  functions  associated  with  the  mRNAs  expressed  in  cells  exposed  to  THz  radiation  and  BH.  The  graph  represents  the  molecular 
and  cellular  functions  relative  to  the  total  number  of  mRNAs  associated  with  each  function.  Twenty-two  biological  functions  are  common  to  THz  and  BH,  four 
are  unique  to  THz  and  one  function  is  unique  to  BH. 


exposure,  respectively,  as  shown  in  Fig.  4(c).  In  both  THz  radia¬ 
tion  and  BH-exposed  cells,  the  number  of  up-regulated  mRNAs 
was  higher  than  that  of  down-regulated  mRNAs.  The  up-regu- 
lated  mRNAs  were  75%  and  85%  at  four  hours  post  THz  radi¬ 
ation  and  BH  exposure,  respectively. 

To  determine  whether  any  of  the  differentially  expressed 
mRNAs  were  scored  in  both  or  solely  in  THz  radiation  or  BH 
exposed  cells,  a  biomarker  comparison  analysis  was  performed 
using  IPA.  This  tool  allowed  us  to  separate  the  identified 
mRNAs  into  three  groups:  unique  to  THz  radiation,  unique  to 
BH,  and  shared  between  THz  radiation  and  BH.  As  shown  in 
the  Venn  diagram  in  Fig.  5(a),  there  were  a  total  of  479  and  620 
mRNAs  that  were  unique  to  THz  radiation  and  BH  exposure, 
respectively.  Additionally,  there  were  52  mRNAs  shared  by 
both  THz  radiation  and  BH  exposures.  The  52  common  mRNAs 
represent  less  than  10%  of  the  combined  total  differentially 
expressed  mRNAs  after  exposure  to  THz  radiation  or  BH. 
Thus,  this  considerable  difference  in  differentially  expressed 
mRNAs  between  THz  and  BH  exposures  suggests  that  the 
intracellular  pathways  triggered  in  either  THz  radiation  or  BH 
might  be  vastly  different. 

The  altered  mRNAs  in  each  group  were  further  separated  by 
their  corresponding  gene  families.  The  pie  charts  in  Fig.  5(b)  il¬ 
lustrate  the  gene  families  related  to  the  unique  and  shared  genes 
as  well  as  the  number  of  genes  for  each  family.  The  gene  fam¬ 
ilies  are  color  coded  as  depicted  in  the  legend.  The  differen¬ 
tially  expressed  genes  belong  to  thirteen  known  families  (cy¬ 
tokines,  enzymes,  G-protein  coupled  receptors,  growth  factors, 
ion  channels,  kinases,  ligand-dependent  nuclear  receptors,  pep¬ 
tidases,  phosphatases,  transcription  regulators,  translation  regu¬ 


lators,  transmembrane  receptors,  and  transporters)  in  addition  to 
genes  that  were  grouped  by  IPA  in  “other”  gene  families.  There 
were  a  few  shared  genes  in  most  gene  families  except  for  cy¬ 
tokines,  peptidases,  growth  factors,  and  translation  regulators 
families,  as  represented  by  the  protruded  sectors  in  Fig.  5(b). 
In  these  four  gene  families,  the  differentially  expressed  genes 
identified  were  all  unique  to  either  THz  radiation  or  BH  gene 
group. 

B.  Biological  Functions  of  the  THz  and  BH  Differentially 
Expressed  mRNAs 

After  identification  of  the  mRNAs  that  were  differentially  ex¬ 
pressed  in  response  to  THz  radiation  and/or  BH,  a  characteri¬ 
zation  of  each  gene’s  primary  cellular  and  molecular  functions 
was  also  performed.  The  bar  graph  in  Fig.  6  represents  the  total 
number  of  genes  associated  with  each  biological  function.  The 
differentially  expressed  genes  in  THz  radiation  and  BH  shared 
22  cellular  and  molecular  functions.  Four  biological  functions 
(energy  production,  free  radial  scavenging,  RNA  damage  and 
repair,  and  vitamin  and  mineral  metabolism)  were  unique  to 
genes  from  the  THz  radiation  exposure,  and  one  function  (pro¬ 
tein  synthesis)  was  unique  to  genes  from  the  BH  exposure. 

C.  Pathway  Analysis  of  Microarray  mRNAs  Data 

IPA  was  also  used  to  determine  the  canonical  pathways  that 
were  affected  in  cells  in  response  to  THz  radiation  and  BH.  54 
metabolic  pathways  and  211  signaling  pathways  were  scored  in 
THz  radiation  exposure,  whereas  45  metabolic  pathways  and 
155  signaling  pathways  were  scored  in  BH  exposure.  After  set¬ 
ting  the  significant  likelihood  cutoff  of  all  canonical  pathways  at 
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Fig.  7.  THz  radiation  and  BH  induce  distinct  sets  of  metabolic  pathways, 
(a)  Significant  metabolic  pathways  and  related  differentially  expressed  genes 
from  THz  radiation  exposure,  (b)  Significant  metabolic  pathways  and  related 
differentially  expressed  genes  from  BH  exposure.  The  metabolic  pathways 
significance  cutoff  was  set  at  threshold  value  (p  —  value  <  0.05).  The 
-log  (p-value)  is  the  maximum  p-value  for  all  of  the  genes  in  each  particular 
canonical  pathway  family.  The  ratio  plot  indicates  the  number  of  the  genes 
expressed  at  each  time  point  relative  to  the  total  number  of  genes  in  that 
particular  metabolic  pathway  group.  Bar  graphs  that  represent  the  magnitude 
of  the  log2  fold-change  in  expression  of  the  genes  associated  with  the  scored 
metabolic  pathways  in  THz  radiation  or  BH  exposure  compared  to  sham  are 
shown  in  (a)  and  (b)  lower  panels. 


the  threshold  value  (p  —  value  <  0.05),  12  metabolic  pathways 
and  1 8  signaling  pathways  were  scored  significant  in  the  THz 
exposure,  as  seen  in  Fig.  7(a)  and  Fig.  8(a),  respectively.  For 


the  BH  exposure,  three  metabolic  pathways  and  13  signaling 
pathways  were  significant,  as  shown  in  Fig.  7(b)  and  Fig.  9(a), 
respectively.  All  the  significantly  affected  canonical  pathways 
in  THz  radiation  and  BH  exposures  were  unique  to  either  expo¬ 
sure.  In  fact,  no  overlap  in  the  metabolic  or  signaling  pathways 
between  THz  radiation  and  BH  exposures  was  observed. 

Next,  the  significant  canonical  pathways  were  sorted  by 
their  respective  categories  in  order  to  get  a  clear  picture  of 
their  roles.  The  metabolic  pathways  associated  with  THz  ex¬ 
posure  are  involved  in  processes  of  degradation,  utilization, 
and  detoxification,  whereas  BH  metabolic  pathways  are  linked 
to  biosynthesis  and  transport.  For  the  THz  radiation  and  BH 
significantly  scored  signaling  pathways,  although  the  signaling 
pathways  were  specific  to  the  type  of  exposure,  most  belong  to 
four  pathway  categories  that  are  present  in  both  THz  radiation 
and  BH.  These  pathway  categories  are  cellular  growth  prolif¬ 
eration  and  development,  intracellular  and  second  messenger, 
organismal  growth  and  development,  cellular  immune  re¬ 
sponse,  and  nervous  system,  as  shown  in  Figs.  8(a)  and  9(a),  top 
panels.  The  signaling  pathway  categories  that  were  exclusive 
to  either  THz  radiation  or  BH  exposure  included  growth  factor, 
cell  cycle  regulation  and  disease-specific  signaling  pathways 
in  THz  radiation  exposure,  and  cardiovascular,  cytokine,  and 
nuclear  receptor  signaling  pathways  in  BH  exposure. 

The  THz  radiation  and  BH  differentially  expressed  genes  that 
are  associated  with  the  significantly  scored  pathways  were  then 
determined.  Twelve  and  seven  differentially  expressed  genes 
were  identified  for  the  metabolic  pathways  of  THz  radiation 
and  BH,  respectively.  Bar  graphs  of  the  effects  of  THz  radi¬ 
ation  or  BH  on  the  expression  profile  of  the  transcripts  from 
these  genes  are  shown  in  the  lower  panels  of  Fig.  7(a)  and  (b). 
These  bar  graphs  represent  the  log2  fold-change  in  expression  in 
THz  radiation  or  BH  exposure  compared  to  sham.  On  the  other 
hand,  46  and  73  differentially  expressed  transcripts  were  identi¬ 
fied  for  THz  radiation  and  BH  signaling  pathways,  respectively. 
Figs.  8(b)  and  9(b)  show  bar  graphs  representing  the  effect  of 
THz  radiation  or  BH  on  the  expression  profile  of  these  mRNAs 
compared  to  sham.  Notably,  the  THz  radiation  and  BH  altered 
genes  did  not  display  a  cellular  component  specific  expression 
but,  rather,  they  included  genes  that  are  expressed  from  all  cel¬ 
lular  locations,  namely,  cytoplasm,  nucleus,  plasma  membrane, 
extracellular  space  and  others  (i.e.,  signalosome,  mitochondrial 
matrix  or  unknown),  as  shown  in  Figs.  8(b)  and  9(b). 

Furthermore,  using  IPA  and  available  web-based  resources, 
the  biological  process  associated  with  the  signaling  path¬ 
ways  related  genes  were  characterized.  Some  of  these  genes 
are  involved  in  precise  biological  processes;  however,  many 
are  linked  to  multiple  key  biological  processes,  as  seen  in 
Figs.  8(b)  and  9(b).  Therefore,  for  uniformity,  these  genes  were 
classified  based  on  top  key  biological  processes  and  overlaps. 
The  results  show  that  THz  radiation  and  BH  exposures  shared 
biological  processes  linked  to  transcription  regulation,  trans¬ 
port,  and  cell  signaling  and  signal  transduction.  Nonetheless, 
the  genes  associated  with  these  processes  were  entirely  different 
between  THz  radiation  and  BH,  except  for  two  genes  (marked 
with  asterisks):  the  transcription  factor  JUN  (Jun  Proto-Onco¬ 
gene)  that  was  induced  10-fold  by  THz  radiation  and  74-fold 
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Fig.  8.  THz  radiation  signaling  pathways  and  related  genes,  (a)  Significant  signaling  pathways  scored  in  THz  radiation  exposure.  The  signaling  pathways  signifi¬ 
cant  likelihood  cutoff  was  set  at  threshold  value  (p  —  value  <  0.05) .  The  —  log  (p-value)  is  the  maximum  p-value  for  all  of  the  genes  in  each  particular  canonical 
pathway  family.  The  ratio  plot  indicates  the  number  of  the  genes  expressed  at  each  time  point  relative  to  the  total  number  of  genes  in  that  particular  signaling 
pathway  group,  (b)  The  differentially  expressed  genes  associated  with  the  THz  significant  signaling  pathways.  Bar  graphs  represent  the  magnitude  of  the  log2 
fold-change  in  expression  of  the  genes  compared  to  sham.  The  THz  radiation  differentially  expressed  genes  were  grouped  according  to  their  top  or  overlapping 
biological  processes  as  shown. 


by  BH  and  the  kinase  NTRK2  (neurotrophic  tyrosine  kinase 
receptor  type  2)  that  was  induced  12-fold  by  THz  irradiation 
and  8-fold  by  BH. 

The  results  also  show  that  the  major  difference  between  ex¬ 
perimental  groups  was  that  THz  radiation  exposure  specifically 
altered  genes  involved  in  biological  processes  of  immune  re¬ 
sponse,  cell  shape  and  adhesion,  and  cytoskeleton  organization. 
Examples  of  the  genes  [shown  in  Figs.  4(a)  and  Fig.  8(b)]  in¬ 
cluded  the  immune  response  genes  [i.e.,  IF4  (interleukin  4)  and 


IGHG1  (immunoglobulin  heavy  constant  gamma  1)],  the  actin 
genes  ACTA1  (alpha-actin-1)  and  ACTA2  (alpha-actin-2),  and 
the  cell  adhesion  genes  ITGB1  (integrin  Beta-1),  ITGB3  (inte¬ 
grin  Beta-3),  CDH19  (cadherin-19),  MYH4  (myosin-4),  RHOB 
(ras  homolog  family  member  B)  and  TFN2  (talin  2),  and  the 
actin  cytoskeleton  organization  molecules  such  as  RHOJ  (ras 
homolog  family  member  J),  FGF  (fibroblast  growth  factor)  7, 
FGF10,  and  WASF  (Wiskott-Aldrich  syndrome-like)).  Alter¬ 
ation  in  the  expression  of  such  genes  underscores  the  signif- 
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Fig.  9.  BH  signaling  pathways  and  related  genes,  (a)  Significant  signaling  pathways  scored  in  BH  exposure.  The  signaling  pathways  significant  likelihood  cutoff 
was  set  at  threshold  value  (p  —  value  <  0.05).  The  —  log  (p-value)  is  the  maximum  p-value  for  all  of  the  genes  in  each  particular  canonical  pathway  family.  The 
ratio  plot  indicates  the  number  of  the  genes  expressed  at  each  time  point  relative  to  the  total  number  of  genes  in  that  particular  signaling  pathway  group,  (b)  The 
differentially  expressed  genes  associated  the  BH  significant  signaling  pathways.  Bar  graphs  represent  the  magnitude  of  the  log2  fold-change  in  expression  of  the 
genes  compared  to  sham.  The  BH  differentially  expressed  genes  were  grouped  according  to  their  top  or  overlapping  biological  processes  as  shown. 


icance  of  signaling  pathways  such  as  the  Actin  Cytoskeleton 
Signaling  and  the  Actin  Nucleation  by  ARP-WASP  complex 
signaling  in  response  to  the  THz  radiation  exposure. 

BH  exposure  [shown  in  Figs.  4(a)  and  9(b)]  specifically  af¬ 
fected  the  expression  of  genes  implicated  in:  1)  response  to  heat 
and  unfolded  protein,  which  included  the  heat  shock  protein 
(HSP)  70  genes  (e.g.,  HSPA1A,  HSPA6,  HSPA1L  and  HSPA2), 
HSP40  genes  (e.g.,  DNAJB  1,  4  and  7),  as  well  as  smaller  HSPs 
like  HSP27  gene  (HSPB1)  and  HSP60  gene  (HSPD1);  2)  prote¬ 
olysis  and  extracellular  matrix  organization  genes  such  as  genes 
from  the  Matrix  Metallopeptidase  (MMP)  family  (e.g.,  MMP8, 
MMP11,  MMP21)  and  from  the  A  Disintegrin  and  Metallopro- 


tease  family  (e.g.,  ADAM  12  and  30);  3)  cell  growth,  prolifer¬ 
ation  or  apoptosis  (e.g.,  DUSP  (dual  specificity  phosphatase)  1, 
10  and  16,  JAK  (Janus  kinase)  3,  SGK  (serum/glucocorticoid 
regulated  kinase)  1 ,  GADD  (growth  arrest  and  DNA-damage- 
inducible)  45B  and  45 G);  4)  transcription  regulatory  molecules 
and  immediate  early  genes,  which  included  JUN,  JUND  (jun  D 
proto-oncogene),  FOS  (FBJ  murine  osteosarcoma  viral  onco¬ 
gene  homolog),  FOSB,  EGR  (early  growth  response)  1,  2,  3 
and  4;  and  5)  nucleosome  assembly  (e.g.,  HIST  (histone  cluster) 
1H1  and  3H3)  and  a  few  more  linked  to  metabolism  of  alcohol, 
carbohydrate  or  lipid,  or  small  molecules.  As  seen  in  Fig.  4(a), 
of  all  the  BH  signaling  pathways  related  genes,  only  HSPA1A, 
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Fig.  10.  THz  radiation  and  BH  alter  the  expression  of  several  miRNAs.  Volcano  plots  of  the  global  expression  profile  of  miRNAs.  The  signal  difference  (log2) 
ratio  is  plotted  versus  the  level  of  statistical  significance  (p-value)  of  all  miRNA  in  the  microarray  for  THz  (upper  panel)  and  BH  (lower  panel)  in  exposed  cells 
compared  to  sham.  Up-regulated  and  down-regulated  miRNAs  with  p  —  values  <0.05  are  depicted  with  green  and  red  triangle  symbols,  respectively.  Insignificant 
genes  are  denoted  with  light  gray  circle  symbols.  A  total  of  107  (47  up-regulated  and  60  down-regulated)  and  90  miRNAs  (33  up-regulated  and  57  down-regulated) 
have  significant  p-values  for  THz  radiation  and  BH  exposure,  respectively. 


HSPA6,  EGR1,  EGR4,  JUN  and  NTRK2  were  also  significantly 
altered  by  THz  exposure. 

D.  THz  Radiation  and  BH  Induce  Differential  Expression  of 
Several  miRNAs 

We  previously  demonstrated  that  dermal  fibroblasts  differ¬ 
entially  express  a  specific  group  of  miRNAs  when  exposed  to 
hyperthermic  stress  [66].  Since  2.52  THz  induces  a  tempera¬ 
ture  increase  in  exposed  Jurkat  cells,  we  sought  to  determine 
whether  THz  radiation  exposure  would  impact  the  expression  of 
miRNAs  similarly  to  thermally-matched  BH  exposure.  miRNA 
microarrays  were  performed  to  identify  differentially  expressed 
miRNAs  in  cells  four  hours  post-exposure.  A  global  view  of 
the  miRNA  profiles  is  provided  as  a  volcano  plot  in  Fig.  10.  In 
these  plots,  log2  signal  difference  ratio  relative  to  the  p-value 
of  all  miRNAs  in  the  microarray  was  mapped  for  THz  (upper 
panel)  and  BH  (lower  panel)  exposures  compared  to  sham. 
Up-regulated  miRNAs  appear  towards  the  top  of  the  plot  (de¬ 
noted  with  green  triangles),  whereas  down-regulated  miRNAs 
appear  towards  the  bottom  of  the  graph  (denoted  with  red 
triangles).  The  significance  value  was  set  at  a  p  —  value  <  0.05. 
Thus,  the  gray  color  circles  represent  the  miRNAs  that  were 
statistically  unchanged.  In  THz  radiation-exposed  cells,  the 
expression  of  1 07  miRNAs  was  significantly  changed,  of  which 
47  were  up-regulated  and  60  were  down-regulated.  Compara¬ 
tively,  in  BH-exposed  cells,  90  miRNAs  displayed  a  significant 
change  in  their  expression  pattern,  with  33  up-regulated  and  57 
down-regulated. 


An  additional  significance  cutoff  was  set  for  the  log2 
signal  difference  ratio,  followed  by  IPA  to  perform  miRNA 
biomarker  comparison  analysis.  With  this  extra  criterion,  the 
expression  of  66  and  53  miRNAs  was  significantly  altered 
(fold  change  >  1.5)  by  THz  radiation  and  BH  exposure, 
respectively.  As  presented  in  the  Venn  diagram  [Fig.  11(a)],  of 
these  miRNAs,  57  were  unique  to  THz  radiation  exposure,  44 
were  unique  to  BH  exposure,  and  nine  were  common  to  both 
exposures.  These  numbers  indicate  that  less  than  17%  of  the 
differentially  expressed  miRNAs  in  THz  radiation  and  BH  were 
shared.  Fold  changes  in  expression  of  these  shared  and  unique 
miRNAs  compared  with  sham  are  shown  in  Fig.  1 1(b)— (d). 
Five  of  the  shared  miRNAs  were  down-regulated  in  both 
THz  radiation  and  BH  exposed  cells,  whereas  four  miRNAs 
(miR-576-5p,  miR-526a,  miR-362-5p,  and  miR-208a-3p)  dis¬ 
played  a  differential  change  in  expression  pattern  between  THz 
radiation  and  BH  exposures.  These  miRNAs  (marked  with 
asterisks)  were  down-regulated  by  THz  radiation  while  induced 
by  BH. 

Next,  to  determine  target  mRNAs  for  the  altered  miRNAs, 
IPA  was  used  to  create  a  miRNA-mRNA  pairing  analysis  for 
the  significantly  differentially  expressed  miRNAs  and  mRNAs 
in  THz  radiation  or  BH  exposure.  From  the  THz-irradiation 
miRNA-mRNA  pairing,  48  of  the  66  altered  miRNAs  have  345 
putative  (predicted  or  experimentally  observed)  target  mRNAs. 
In  the  BH  exposure  miRNA-mRNA  pairing,  there  were  42  of 
the  53  altered  miRNAs  that  have  438  putative  target  mRNAs. 
As  expected,  a  number  of  these  miRNAs  have  multiple  sets 
of  target  mRNAs.  As  an  example,  Table  I  presents  the  result 
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Fig.  1 1 .  Unique  and  shared  miRNAs.  (a)  Venn  diagram  indicating  the  number  of  unique  and  shared  miRNAs  exhibiting  a  significant  level  of  differential  expression 
in  THz  radiation  and  BH  exposures.  The  significant  cutoffs  were  set  as  fold  change  value  >1.5  and  a  p  —  value  <  0.05.  (b)  Bar  graph  representing  the  fold 
change  in  expression  of  the  shared  miRNAs  in  THz  radiation  and  BH  exposure  compared  to  sham,  (c)  Bar  graph  representing  the  fold-change  in  expression  of  the 
unique  miRNAs  in  THz  radiation  compared  to  sham,  (d)  Bar  graph  representing  the  fold-change  in  expression  of  the  unique  miRNAs  in  BH  exposure  compared 
to  sham. 


from  a  pairing  of  each  of  the  nine  THz  and  BH-shared  miRNAs 
[Fig.  1 1  (b)]to  all  of  the  significantly  differentially  expressed 
mRNAs  in  THz  or  BH  exposure.  In  this  table,  THz  putative 
target  mRNAs  that  were  unique  to  THz  or  BH  exposure  as  well 
as  those  that  were  shared  in  both  exposures  are  highlighted. 
As  shown,  six  of  these  miRNAs  (miR-1827,  miR-208a-3p, 


miR-296-5p,  miR-3617-5p,  miR-517a-3p,  and  miR-576-5p) 
have  one  or  few  common  target  mRNAs.  Of  the  remaining 
miRNAs,  miR-362-5p  and  miR-526a  have  a  totally  different 
set  of  target  mRNAs  in  THz  radiation  versus  BH  exposure 
whereas  miRNA  let-7a-5p  has  putative  targets  only  in  the  THz 
radiation  differentially  expressed  mRNAs. 
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TABLE  I 

Comparison  of  THz  and  BH  Target  mRNAs  for  Each  Shared  miRNA.  Each  THz  Radiation  and  BH-Shared  miRNAs  was  Paired  With  the 
Discovered  Differentially  Expressed  mRNAs  From  THz  or  BH  Exposure  Using  IPA.  (-)  Means  no  Putative  Target  mRNAs  Were  Found 


miRNA 

THz  Putative 

Target  mRNAs 

THz  &  BH 
Target 
mRNAs 

BH  Putative 

Target  mRNAs 

let-7a-5p 

ACTA1,  AG04,  BTN2A1,  COL1  Al,  CPM, 
CYP19A1,  DCLK3,  DSP,  DUSP16,  DZIP1, 
EPHA3,  FBX032,  FZD4,  GABBR2,  GHR, 
HAS2,  HMGA2,  ITGB3,  LIMD1, 

LINGO  1,  MSR1,  NOS1,  OTUB2, 

PCYT1B,  PLD5,  PLXNA4,  RHOB, 
SEMA4G,  SFT2D3,  SGCD,  SLC4A7, 
SMUG1,  SULF2,  SYT11,  TTLL4,  WASL, 
ZBTB10,  ZNF493 

miR-1827 

ADAMTS4,  ARSB,  ARSD,  BLOC  1  S3, 
CACNB2,  CALC  A,  CCDC68,  CD300LF, 
CNN3,  CORIN,  COX6B2,  CPM, 

CYP19A1,  DEPDC1,  ECE1,  FAM53B, 
FCH02,  HOXC9,  KHK,  KLHL1,  KY, 
LAMC2,  LANCL2,  LPCAT2,  LRRC25, 
LRRC7,  MAP3K15,  MDH1B,  NOS1, 

NPR3,  PAPPA2,  PLEKHA6,  PLXNA4, 
PRRG1,  PVRL4,  PYCRL,  RAPGEF5, 
SAMHD1,  SIRT5,  SLC10A1,  SMIM10, 
SPATA18,  SUSD5,  SYT9,  TAT, 

TNFRSF9,  TRAK1 

UBE3B 

ADRA1A,  AHNAK,  APOL6,  CACNA1G,  CBFA2T2, 
CNFN,  CYB561,  DAND5,  ESR1,  FAM160B2,  FAM217B, 
FAM46A,  FGF18,  FGFRL1,  GAL3ST1,  HEY1,  IGF1, 
INPP5B,  JAK3,  JPH3,  KIAA1755,  MUL1,  MUM1,  NRG1, 
PIGV,  PRKAA2,  PSME4,  RBPMS,  RHO,  SCML4, 
SERPLNH1,  SGPL1,  SLC30A2,  SMCR7,  SMIM6,  SOCS3, 
SPPL3,  STRAP,  TAD  A3,  THBS1,  THEM4,  TMC5, 
TMED10,  TMPRSS13,  TRA2B,  VIPR2,  XG,  ZBTB3, 
ZC3H12A,  ZCCHC4,  ZNF557,  ZNF585B,  ZNF689, 
ADAM30 

miR-208a-3p 

EPHA3,  EPHA5,  FZD4,  SOS2,  ARSD 

BHLHE41 

CD14,  CELSR1,  CNKSR2,  COL4A3,  EYA4, 

HIST2H2BE,  MLF1,  MTMR3,  NAV2,  PHLDA2,  PURB, 
SRCIN1,  YTHDC1,  ZEB2 

miR-296-5p 

ABCF2,  ACSL5,  ADRA1D,  AGFG2, 

CD44 

APOL6,  ATP2B3,  AZU1,  Cllorf87,  C19orf26,  CAD, 

ARHGEF4,  ARSB,  BLOC  1  S3,  C9orfl73, 

GPR173 

CADM3,  CELSR1,  CERS1,  CHAC1,  CMBL,  DNAJB1, 

CD300LF,  CNN1,  COX18,  ECE1,  EHMT1, 

HES5 

EPN3,  FAM160B2,  FAM169B,  FAM19A5,  FBN2, 

F9,  FAM53B,  FBX031,  FBX032,  FCRL5, 

MEOX1 

FCAMR,  FGFBP2,  FGFRL1,  FMNL3,  FZD10, 

FGF7,  FGFR3,  FZD4,  GABBR2,  GLI2, 
H1FX,  IFNK,  KHK,  KIF1A,  LINGOl, 
LITAF,  LRCH4,  MCAM,  MRC2, 

NBEAL2,  NCAM1,  NOS1,  OTUB2, 

PIEZ02,  PYCRL,  RAB1  IB,  SLRT5, 

SLC1A7,  SMOC2,  SRRM3,  SSTR1, 

SYDE2,  SYT11,  TC2N,  TNFAIP2,  TTLL4, 
VWA1,  ZNF555 

NUMA1 

GADD45G,  GAL3ST1,  GIGYF1,  HELZ2,  IL11,  IL27, 
JPH3,  KIAA2018,  KLF10,  KLHDC10,  LRRC10B, 
LSMEM1,  LTB4R,  MAP3K8,  MAP7D2,  MEGF6, 

MMP1 1,  NR2F6,  NR4A3,  PHOSPHOl,  PNPLA7, 

POU3F1,  RGS11,  RHBDL3,  SCAMP5,  SNAI1,  SPPL3, 
SRCIN1,  ST3GAL3,  TBC1D13,  TECR,  TLN1, 
TMEM132B,  TMEM158,  TPPP,  USP22,  VAC  14,  WDR20, 
XG,  ZBTB42,  ZNF585B,  ACCS 

miR-3617-5p 

EPHA3,  MFAP5,  NPL,  PGR,  RAPGEF6, 

NPAS3 

CDH8,  EIF4A2,  GLCCI1,  HOXB2,  HPCAL4,  HTR2A, 

RNF17,  SYDE2,  TAC1,  TRAK1,  ZNF677, 
Clorfl73 

CD44 

MS4A6A,  NEFH,  PRKAA2,  SLC22A16,  YTHDC1 

miR-362-5p 

Clorfl73,  C7,  CYLD,  DAZ2,  DAZL, 

NPR3,  OLIG3,  RAPGEF6,  RIPK2, 

SLITRK6,  TRAK1,  C18orf54 

C15orf41,  CLCA4,  DLL1,  DNAJB1,  GPR98,  MAP3K8, 
MLF1,  MRC1,  SLC35G3 

miR-517a-3p 

Clorfl73,  HPGDS,  MYH4,  RALGAPA1, 
SMIM10 

BHLHE41 

ANKRD29,  CDKN1C,  H2BFM,  IGF1,  IRF2BP2, 
LGALS13,  MS4A4A,  NPAS4,  PDE4D,  SIGLEC5,  TNKS, 
WDR91 

miR-526a 

CASD1,  CREBRF,  CYLD,  ELL2,  FGF7, 
JPH4,  MCF2L,  MFAP5,  MSR1,  NPL, 
PRRG1,  RALGAPA1,  SFT2D3,  SYDE2 

ATF3,  CCNL1,  CDH8,  CXCL1 1,  DOCK5,  EFNB2, 
MORC1,  MRC1,  SLC1A2,  TNFRSF19,  TNFSF10 

miR-576-5p 

CD36,  FCH02,  HOPX 

BHMT 

Cllorf87,  DNAJB7,  GJA5,  ITPR1,  KLF10,  MGARP, 
NAMPT 

Each  THz  radiation  and  BH-shared  miRNAs  was  paired  with  the  differentially  expressed  mRNAs  from  THz  or  BH  exposure  using  IPA. 
(-)  means  no  putative  target  mRNAs  were  found. 


IV.  Discussion  and  Conclusion 

In  the  present  study,  we  aimed  to  perform  a  full  transcrip¬ 
tional  response  of  mRNAs  and  miRNAs  and  used  these  data  to 
identify  the  intracellular  canonical  pathways  that  are  triggered 
in  THz-irradiated  cells.  We  specifically  sought  to  answer  the  fol¬ 
lowing  fundamental  question:  Do  cells  exposed  to  high-power 
2.52  THz  radiation  exhibit  similar  or  different  specific  cellular 
and  molecular  events  compared  to  cells  exposed  to  an  equiv¬ 
alent  dose  (temperature  rise  and  exposure  time)  conventional 
thermal  BH? 


In  our  experiments,  a  40  min  exposure  of  Jurkat  cells  to  a  CW, 
single-frequency,  2.52  THz  produced  about  a  6  °C  rise  in  tem¬ 
perature.  Therefore,  we  included  a  matched  44  °C  BH  control 
for  comparison.  It  is  interesting  to  note  that  our  results  showed 
that  cells  exposed  to  2.52  THz  radiation  differentially  expressed 
a  specific  set  of  mRNAs  and  miRNAs  than  those  differentially 
expressed  in  BH-exposed  cells.  In  fact,  less  than  10%  and  17% 
of  the  altered  mRNAs  and  miRNAs,  respectively,  were  shared 
by  THz  and  BH  exposures.  As  a  result,  we  observed  that  these 
differentially  expressed  genes  involve  specific  metabolic  and 
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signaling  pathways  in  THz  radiation  versus  BH  exposure  when 
compared  to  the  sham  control. 

It  is  known  that  a  temperature  increase  due  to  heating  can 
cause  intracellular  and  membrane  proteins  to  unfold,  entangle, 
and  aggregate  [46]— [53].  Such  modifications  can  preclude  cells 
from  functioning  properly.  Therefore,  in  order  to  overcome 
such  stress  and  ensuing  pleotropic  effects,  cells  utilize  their 
cellular  defense  and  adaptation  mechanism,  commonly  referred 
to  as  the  cellular  stress  response  (CSR).  The  CSR  is  rapidly 
activated  in  response  to  stressors  such  as  heat,  and  a  hallmark 
of  this  response  is  a  change  in  expression  of  hundreds  of  genes 
[46].  The  stress-inducible  genes  encode  for  a  multitude  of  pro¬ 
teins  involved  in  many  important  functions  including:  cellular 
survival,  cell  cycle  control,  and  apoptosis;  molecular  chaper¬ 
ones,  with  the  predominant  class  being  heat  shock  proteins 
(Hsps);  regulatory  proteins,  including  transcription  regulatory 
molecules,  cytokines,  and  immediate  early  genes  (IEGs); 
cell  organization  and  structure,  such  as  proteins  required  to 
sustain  actin  cytoskeleton;  proteolysis  (i.e.,  proteases);  DNA 
sensing  and  repair;  metabolism;  and  membrane  regulation  and 
transport  [67],  [68].  In  this  regard,  it  is  not  surprising  that 
BH  elicited  a  heat  stress  response  in  the  exposed  cells,  which 
was  manifested  through  a  significant  increase  in  the  mRNAs 
encoding  for  molecular  chaperoning  proteins  including  the 
well  characterized  HSP70  (HSPA1A,  HSPA6,  HSPA1L  and 
HSPA2)  and  HSP40  (DNAJB  1,  4  and  7),  as  well  as  smaller 
proteins  like  HSPB1  and  HSPD1.  Contrarily,  although  THz 
radiation  induced  a  marked  temperature  rise  (in  the  order  of  few 
degrees),  it  only  caused  an  induction  in  HSPA1A  and  HSPA6. 
The  expression  of  all  the  other  mRNAs  encoding  for  heat 
stress  and  chaperoning  were  not  significantly  changed  by  THz 
radiation,  suggesting  that,  in  our  condition,  2. 2 5 -THz  radiation 
does  not  appear  to  be  directly  damaging  intracellular  proteins. 

Furthermore,  BH  exposure  altered  the  expression  of  several 
other  CSR  gene  functions  such  as  proteolysis  (MMP8,  MMP1 1 , 
MMP21,  ADAM  12  and  30),  extracellular  matrix  organization, 
cell  growth,  proliferation,  or  apoptosis  (e.g.,  DUSPs  JAK3, 
SGK1,  GADD45B  and  45G)  in  addition  to  transcription  regu¬ 
latory  molecules  and  immediate  early  genes.  THz  radiation,  on 
the  contrary,  had  no  significant  effect  on  the  expression  of  most 
of  these  genes.  THz  radiation  exposure  altered  genes  of  the  im¬ 
mune  response  (IL4  and  IGHG1)  and  the  expression  of  several 
genes  involved  in  cell  shape,  adhesion,  and  cytoskeleton  organ¬ 
ization  (e.g.,  ACTA1  and  ACTA2,  ITGB1,  ITGB3  CDH19  and 
RHOB),  which  were  not  significantly  changed  in  BH. 

Our  data,  show  that  there  is  a  discrepancy  between  the  genes, 
metabolic  and  signaling  pathways  enriched  for  THz  radiation 
versus  BH  exposure.  They  also  provide  evidence  that  THz  ra¬ 
diation  did  not  alter  the  expression  of  major  heat  shock  chaper¬ 
oning  and  proteolysis  proteins  that  were  affected  in  BH.  There¬ 
fore,  the  results  imply  that  2. 25 -THz  radiation  does  not  appear  to 
directly  denature  or  significantly  damage  intracellular  proteins 
and  suggest  that  the  heating  regimen  induced  by  THz  radiation 
is  different  from  that  of  the  conventional  thermal  BH.  More¬ 
over,  since  the  dynamics  and  vibrational  motion  of  interfacial 
water  differ  considerably  from  that  of  bulk  water,  and,  given 
that  interfacial  water  absorption  is  maximal  at  ^2.4-2. 5  THz 
[12],  we  speculate  that  electromagnetic  radiation  at  a  frequency 


of  2.5  THz  is  preferentially  coupled  to  and  selectively  absorbed 
by  interfacial  water  molecules  in  living  cells.  This  selective 
coupling  mechanism  might  cause  appreciable  local  tempera¬ 
ture  rises  on  the  surface  of  biomolecules,  triggering  specific  cel¬ 
lular  responses  not  observed  with  the  BH  protocol.  Nonetheless, 
using  a  different  heating  protocol  than  the  conventional  thermal 
BH,  such  as  using  a  near  IR  beam,  to  generate  a  same  tempera¬ 
ture  rise  as  THz,  would  be  needed  to  validate  our  speculations. 

Inherent  gene  regulation  is  fundamental  to  the  versatility 
and  adaptability  of  cells.  Therefore,  a  change  in  such  processes 
by  an  external  stimulus  such  as  THz  radiation  could  lead  to 
significant  biochemical  or  physiological  implications.  Hence, 
the  use  of  THz  radiation  as  a  non-contact  tool  for  the  selective 
targeting  of  specific  gene  programs  would  have  a  considerable 
impact  in  biotechnology  and/or  medicine.  In  such  a  scenario, 
THz  radiation  could  potentially  be  used  to  selectively  regulate 
key  genes  involved  in  specific  pathways  and  cellular  processes. 
Few  examples  include  cellular  reprogramming,  cellular  pre¬ 
conditioning,  and  nerve  stimulation  or  repression  [38],  [39], 
[69],  [70].  A  recent  study  showed  that  exposure  of  mouse  stem 
cells  to  2. 25 -THz  radiation  caused  gene  specific  transcriptional 
alterations  and  accelerated  stem  cell  differentiation  toward 
adipocytes  [38],  [39].  Furthermore,  another  study,  showed  the 
stimulating  effect  of  THz  radiation  on  nerve  cells  [70]. 

Our  study  provides  valuable  new  insights  that  present  a 
much  clearer  picture  of  intracellular  canonical  pathways  that 
are  specifically  triggered  in  human  cells  exposed  to  high-power 
CW  2.52  THz.  Further  investigations  involving  other  cell  lines 
and  primary  cells  with  a  variation  in  the  THz  exposure  regimes, 
time  courses,  size,  shape  and  doses  of  the  THz  beam  and  studies 
using  a  wide  range  of  THz  radiation  frequencies  are  needed  to 
validate  our  conclusions. 
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